Exact results for Casimir interactions between dielectric bodies: The weak-coupling or 

van der Waals Limit 
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In earlier papers we have applied multiple scattering techniques to calculate Casimir forces due 
to scalar fields between different bodies described by delta function potentials. When the coupling 
to the potentials became weak, closed-form results were obtained. We simplify this weak-coupling 
technique and apply it to the case of tenuous dielectric bodies, in which case the method involves 
the summation of van der Waals (Casimir-Polder) interactions. Once again exact results for finite 
bodies can be obtained. We present closed formulas describing the interaction between spheres and 
between cylinders, and between an infinite plate and a retangular slab of finite size. For such a slab, 
we consider the torque acting on it, and find non-trivial equilibrium points can occur. 

PACS numbers: 03.70.+k, 03.65.Nk, 11.80.La, 42.50.Lc 



The subject of the perturbation of the quantum vac- 
uum energy by material bodies, the so-called Casimir 
effect, has a long history For only a limited num- 
ber of situations, possessing a high degree of symmetry 
(the interaction energy between infinite parallel planes 
0, the self-energy of spheres [|| and cylinders have 
exact (numerical) calculations been possible. The ex- 
periments carried out during the last decade or so have 
typically been for a spherical surface above a plane sur- 
face. (For a review of the experimental situation, see 
Ref. |5[.) Because the forces for that geometry could 
not be calculated exactly, comparision with theory was 
made using the proximity force theorem (PFT), which 
allows one to compute the force between curved surfaces 
which are nearly touching Q . However, there is no well- 
defined "proximity force approximation" that allows one 
to calculate corrections to the PFT in powers of the ratio 
of the separation distance to the radius of curvature of 
the surface. Because the precision of the experiments is 
now approaching 1%, such corrections may become im- 
portant. 

Various interesting developments have improved the 
theoretical situation. For example, there has been no- 
table progress in developing the numerical Monte-Carlo 
worldline method of Gies and Klingmuller The dif- 
ficulty with this technique lies in the statistical limita- 
tions of Monte Carlo methods and in the complexity of 
incorporating electromagnetic boundary conditions. Op- 
tical path approximations have been studied extensively 
for many years, with considerable success Q. However, 
there always remain uncertainties because of unknown 
errors in excluding diffractive effects. Direct numerical 
methods 0|, based on finite-difference engineering tech- 
niques, may have promise, but the requisite precision 
of 3-dimensional calculations may prove challenging. A 
methodology which is, in principle, exact is the multiple 
scattering formalism, which dates back at least into the 



1950 s [lol . 11, 12, lj| . For more complete references see 
Ref. [ljf 

Previous work on this multiple scattering technique, 
which has been brought to a high state of perfection by 
Emig et al. [l3| . has concentrated on numerical results 
for the Casimir forces between conducting and dielectric 
bodies such as spheres and cylinders. Recently, we have 
noticed that the multiple-scattering method can yield ex- 
act, closed- form results for bodies that are weakly cou- 
pled to the quantum field |l4|. This allows an exact as- 
sessment of the range of applicability of the PFT. The 
calculations there, however, as those in recent extensions 
of our methodology fl5| . have been restricted to scalar 
fields with 5-function potentials, so-called semitranspar- 
ent bodies. (These are closely related to plasma shell 
models.) Here we remedy that defect by making the ex- 
tension of the formalism to electromagnetism, and the 
bodies are, correspondingly, characterized by a permit- 
tivity or dielectric constant e. Strong coupling would 
mean a perfect metal, e — > oo, while weak coupling means 
that e is close to unity. 

In this Letter we will briefly review the formalism, and 
show that it is precisely equivalent to summing Casimir- 
Polder or van der Waals (vdW) potentials. Exact results 
have been found in the past in such summations, for ex- 
ample for the self-energy of a dilute dielectric sphere [l6| 
or a dilute dielectric cylinder [l7|. Thus it is not sur- 



prising that exact results for the interaction of different 
dilute bodies can be obtained. It is only surprising that 
such results were not found much earlier. (We note that 
the additive approximation has been widely used in the 
past, for example see Ref. (l8| . but here the method is 
exact.) We will consider a slab of finite extent above an 
infinite plane, and calculate the force and torque on the 
slab. With the center of mass fixed, we find generically 
the shortest side of the slab aligns with the plate, but for 
sufficiently square slabs nontrivial equilibrium points can 
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FIG. 1: Slab (with dielectric constant £2) of finite cross section 
A but of infinite extent in the z-direction, a distance a above 
an infinite plate (with dielectric constant ei) extending to 
2 = —00. 



be found. We will also compute the force between spheres 
and parallel cylinders. Since the results are exact, we can 
quantify the deviation from the PFT. 

For electromagnetism, we can start from the formalism 
of Schwinger [19(, which is based on the electric Green's 
dyadic T. Just as in the scalar case, the vacuum energy 
for a static configuration, existing for a time r, is given 
by 
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and again, in precise analogy with the scalar case, the 
interaction in lowest order between two non-overlapping 
potentials V% and V2 is 



E= — Tr^r T/ 2 r . 

AT 
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Here, the trace is over both vector indices and space-time 
coordinates, and T is the free Green's dyadic. 
Now an easy calculation shows that 
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where the bodies, which are presumed to be composed of 
uniform material filling nonoverlapping volumes v\ and 
i>2, respectively, are characterized by dielectric constants 
El and £2, both nearly unity. This is the famous Casimir- 
Polder potential [H. 

Consider first a dilute dielectric slab, a distance a 
above a dilute dielectric plate filling the half space z < 0. 
Let the slab have cross section A, and extend from z = a 
to infinity, as shown in Fig. [TJ Because the plate is of 
infinite extent in the x-y plane, it is easy to carry out the 
integrals in Eq. with the result being precisely the 
dilute limit of the Lifshitz formula: 
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from which the force per area is obtained by taking the 
negative derivative with respect to a. There is no correc- 
tion due to the finite size of the slab. This is presumably 
a consequence of the fact that for weak coupling, multiple 
scattering reduces to the two-scattering approximation. 

Instead of integrating over z for the slab, we could have 
considered a slab of thickness dz and 3X6 el v4_ ; 3 distance 
z above the infinite plate, which has the energy 
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This allows us to consider an arbitrarily shaped body 
above the infinite plate. For example, we can immedi- 
ately find the energy of a dilute sphere of radius a, the 
center of which is a distance Z above the plate, Z > a. 
The energy is 



E = 



23 . , . v 

( £l - l)(e 2 - 1) 



1 



640tt 



z 4 {i-a 2 /z 2 y 



(0) 



where v is the volume of the sphere. When the sphere 
nearly touches the plate, S = Z — a <C a, we recover the 
proximity force theorem: 

It is interesting to compare the correction implied by 
the exact sphere-plate energy ([6]) to the proximity force 
theorem result (JT]) 



(8) 



with v = 1. This is to be contrasted with the results 
found by Wirzba [l(| and by Bordag and Nikolaev 
for a scalar field with Dirichlet boundary conditions, v — 
5/ir 2 -l/3 = 0.173; Maia Neto et al. find for perfectly 
conducting boundary conditions for electromagnetic field 
fluctuations v « 1.4. So our dilute dielectrics have an 
intermediate behavior. 

Now, let us return to the slab geometry, but with finite 
size in all directions, and tilted with respect to the infinite 
plate. For simplicity, we will consider only a tilt 6 in 
the y-z plane. See Fig. [2l The center of the body is a 
distance Z above the plane. The result of integrating 
over the body coordinates is the interaction energy 
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FIG. 2: Rectangular solid (with dielectric constant £2) of side 
a, b, and c (perpendicular to the plane, not shown) a distance 
Z above an infinite plate (with dielectric constant £1) extend- 
ing to z = —00. The side b makes an angle 9 with respect to 
the plate. 



interesting here is the torque exerted about the center of 
mass of the slab, r = —dE/88. Generically, the shorter 
side wants to align with the plate, which is obvious geo- 
metrically, since that (for fixed center of mass position) 
minimizes the energy. However, if the slab has square 
cross section, the equilibrium position occurs when a cor- 
ner is closest to the plate, also obvious geometrically. But 
if the two sides are close enough in length, a nontrivial 
equilibrium position between these extremes can occur. 

The situation is illustrated nicely in Fig. [3] Consider a 
slab of aspect ratio b/a < 1. When one corner of the slab 
just touches the plate, the center of mass is Z = D = 
y/a 2 + b 2 /2 and the equilibrium angle satisfies tan0 = 
b/a. For large separation, Z > Zq, where 



(11) 





FIG. 3: The stable equilibrium angle # cq of a slab above an 
infinite plate for given b/a ratios 0.95, 0.9, and 0.7, respec- 
tively given by solid, dashed, and dot-dashed lines. For large 
enough separation, the shorter side wants to face the plate, 
but for Z < Zo the equilibrium angle increases, until finally 
at Z = D = \/ a 2 + b 2 /2 the slab touches the plate at an angle 
8 = arctan6/a, that is, the center of mass is just above the 
point of contact, about which point there is no torque. 



and v is the volume of the slab. This always represents 
an attractive force between the slab and the plate. Most 



the stable equilibrium configuration occurs when 9 = 0, 
that is, when the shortest face b is closest to the plate. 
For distances between Z and D, the stable equilibrium 
angle is intermediate between and tan6/a < 7r/4, as 
shown in Fig. [3] For a = b the stable equilibrium position 
is always at 8 = 7r/4, and for b = it is at 8 = 0. 

Obviously the method described here can be used to 
examine interactions between bodies of arbitrary shape. 
We will content ourselves here with discussing interac- 
tions between cylinders and spheres. 

First, we remark that it is extremely easy to use this 
method to reproduce the results given in Ref. 14 1 
parallel cylinders described by (5-function potentials 
is a bit more complicated to do the calculation for dielec- 
tric cylinders, of radii a and b, separated by a distance 
R > a + b, because we have to integrate over the radii of 
the cylinders. With the 1/r 7 van der Waals potentials, 
the angular integrals can be carried out explicitly, leav- 
ing us with a fairly complicated function of the radial 
coordinates on each cylinder. To proceed, we expand the 
integrand of those radial integrals in r/R, r'/R, and in- 
tegrate term by term, rdr J Q r'dr'. Once more, it is 
remarkable that we can explicitly carry out the sum: 
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For two spheres, with center-to-center distance R and angular integration leads to elliptic integrals of compli- 
radii a and b respectively, the calculation is a bit more cated argument, so it is difficult to proceed using closed- 
complex than for cylinders. The reason is that even one form expressions. Therefore, we expand right at the be- 
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ginning, before carrying out any integrations. We then 
integrate over polar and azimuthal angles, and then over 
the radial coordinates. We can obtain a power series 
expansion in powers of a/ R and b/R, which it is possi- 
ble to sum, once the power series coefficients are identi- 



fied. Again we test the procedure by verifying that we 
reproduce the results found previously for semitranspar- 
ent spheres 14J . Mathematica is actually able to sum the 
resulting series, with the following result for the energy 



E 



23 ( ei -l)( e2 _l) 



1920tt 



R 



In 



4,2 2>4 , >0 

— a b —a o -j-o 



3a -14aV+3fc 4 

w 



o a 2 +b 2 



\ R ) 



R 2 



i-b \ 
R ) 



(13) 



This expression, which is rather ugly, may be verified to 
yield the proximity force theorem: 
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5 = R — a — b <C a,b. It also, in the limit b — > oo, R 
with R — b = Z held fixed, reduces to the result (j 
the interaction of a sphere with an infinite plate. 

In this Letter we have shown that the general method- 
ology of the multiple scattering formulation for weak cou- 
pling becomes the pairwise summation of van der Waals 
energies between the molecules that make up the dilute 
dielectrics. Such summations have previously been car- 
ried out both for the interaction energy between parallel 
plates [19jand the self-interactions of spheres [16[ and 
cylinders [17] . but it was apparently not recognized that 
it was easy to obtain exact closed-form energies for many 
interesting situations. Since the calculations here refer 
to electromagnetic field fluctuations, causing forces and 
torques on dielectric bodies, these are of far more rele- 
vance than our earlier exact results for weak semitrans- 
parent bodies. These results raise the intriguing possi- 
bility that maybe even in strong coupling, for example, 
for conducting bodies, exact results may be obtainable. 
In any case, the results presented here form the basis for 
a laboratory to study edge effects and other finite-size 
phenomena which have proved elusive in the past. 
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